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Abstract 

■The resistivity in magnetic materials has been theoretically shown to depend on the spin-spin correlation function which 
"in turn depends on the magnetic-field, the density of conduction electron, the magnetic ordering stability, etc. However, 
these theories involved a lot of approximations, so their validity remained to be confirmed. The purpose of this work is to 
show by newly improved extensive Monte Carlo (MC) simulation the resistivity of the spin resistivity from low-T ordered 
phase to high-T paramagnetic phase in ferromagnetic and antiferromagnetic films. We take into account the interaction 
between the itinerant spins and the localized lattice spins as well as the interaction between itinerant spins themselves. 
We show that in ferromagnets the resistivity shows a sharp peak at the magnetic phase transition in agreement with 
previous theories in spite of their numerous approximations. Resistivity in antiferromagnets on the other hand shows no 
peak for the SC, BCC and diamond lattices. Discussion on the origin of these resistivity behaviors is given. 
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1. Introduction 

The magnetic resistivity has been extensively studied 
by both theories and experiments in the last fifty years. 
Experiments have shown that the resistivity indeed de- 
pends on the itinerant spin orientation and the lattice 
spin ordering [ll-[6|. At low temperature (T), the main 
magnetic scattering is due to spin- wave excitations 0, HI . 
^he resistivity is proportional to T^. However at higher 
T the spin-wave theory is not valid, such a calculation of 
,the resistivity is not possible, in particular in the critical 
■region around the Curie temperature Tc in simple ferro- 
magnets, let alone other complicated magnetic orderings. 
Experiments on various magnetic materials have found 
in particular an anomalous behavior of the resistivity at 
the critical temperature where the system undergoes the 
ferromagnetic-paramagnetic phase transition [2-6]. Very 
recent experiments such as those performed on ferromag- 
netic SrRuOs thin films[9], superconducting BaFe2As2 sin- 
gle crystals[10j, Lai _^^MnOq[llj. Mni_^Cr^Te[12] and 
other compounds 1 13l4l6| show different forms of anomaly 
of the magnetic resistivity at the transition temperature, 
de Gennes and Friedel's first explanation in 1958 [17J for 
the resistivity behavior near Tc was based on the inter- 
action between the spins of conduction electrons and the 
lattice spins. The resistivity was thus expected to depend 
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strongly on the spin ordering of the system. They have 
suggested that the magnetic resistivity is proportional to 
the spin-spin correlation, therefore it should behave as 
the magnetic susceptibility with a divergence at Tc due 
to "long-range" fiuctuations of the magnetization. Other 
authors [18-20] subsequently suggested that the shape of 
the resistivity results mainly from "short-range" correla- 
tion at T > Tc. Fisher and Langer [19] have shown in 
particular that the form of the resistivity cusp depends on 
the correlation range. To see more details on the role of the 
spin-spin correlation, we quote a work by Haas [21] and a 
more recent work of Kataoka[22] where the spin-spin corre- 
lation function has been calculated. Recently, Zarand et al 
[23] have used the picture that the itinerant spin is mainly 
scattered by impurities which are characterized by a "lo- 
calization length" in the sense of Anderson's localization. 
They found that the peak's height depends on this localiza- 
tion length. Note that since the giant magnetoresistance 
(GMR) was discovered experimentally twenty years ago in 
magnetic multilayers [ij, [111, intensive investigations on 
the spin resistivity, both experimentally and theoretically, 
have been carried out |26|, |27| . The " spintronics" was born 
with a spectacular rapid development in relation with in- 
dustrial applications. For recent overviews, the reader is 
referred to Refs. [28| and [29|. In spite of these inten- 
sive investigations, except our works [30[ Hll, there have 
been no Monte Carlo (MC) simulations performed regard- 
ing the temperature dependence of the spin transport. In 
these works , we have investigated by MC simulations the 
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effects of magnetic ordering on the spin current in mag- 
netic multilayers. 

In this paper we improve our previous MC simulations 
to study the transport of itinerant electrons in ferromag- 
netic and antiferromagnetic crystals. We use the Ising 
model and take into account interactions between lattice 
spins and itinerant spins. We show that in ferromagnets 
we obtain, with our new MC averaging method, much bet- 
ter results for the magnetic resistivity which shows a huge 
peak at the transition temperature while in antiferromag- 
nets, the resistivity does not show such a peak. 

The paper is organized as follows. Section 2 is de- 
voted to the description of our model and the rules that 
govern its dynamics. In section 3, we describe our MC 
method and discuss the results we obtained for ferromag- 
nets. Results on antiferromagnets are shown in Section 4. 
Concluding remarks are given in Section 5. 

2. Model 

2.1. Interactions 

We consider a thin film of FCC structure with two 
symmetrical (001) surfaces. The total number of cells is 
Nx X Ny X Nz where each cell has four spins. Spins lo- 
calized at FCC lattice sites are called "lattice spins" here- 
after. They interact with each other through the following 
Hamiltonian 

H^ = -J^S,-S„ (1) 

where is the Ising spin at lattice site z, XI (ij) iiidi- 
cates the sum over every nearest-neighbor (NN) spin pair 
(Si,Sj), J being the NN interaction. We consider in this 
paper both J > (ferromagnets) and J < (antiferro- 
magnets). 

We consider a flow of itinerant spins interacting with 
each other and with the lattice spins. The interaction be- 
tween itinerant spins is defined by 



(2) 



where is the itinerant Ising spin at position r^, and ^(^^ 
denotes a sum over every spin pair (s^ , s j ) . The interaction 
Kij depends on the distance between the two spins, i.e. 
Tij = |r — Tj |. A specific form of Kij will be chosen below. 
The interaction between itinerant spins and lattice spins 
is given by 



(id) 



(3) 



where the interaction lij depends on the distance between 
the itinerant spin and the lattice spin S^. For the sake 



of simplicity, we assume the same form for Kij and lij^ 
namely. 



Koex.p{-rij) 
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(4) 
(5) 



where Kn and In are constants. 



2.2. Monte Carlo Method 

Before calculating the resistivity, we determine the crit- 
ical temperature Tc below which the system is in the or- 
dered phase using Eq. ([T]). To this end, we perform 
standard Metropolis MC simulations to determine various 
physical quantities at different T[32]. 

Once the lattice is equilibrated at a given T, we inject 
No itinerant spins into the system. The itinerant spins 
move in the x direction under the effect of an electric field. 
We use the periodic boundary conditions to ensure that the 
average density of itinerant spins remains constant with 
evolving time (stationary regime). 

Note that unlike in the previous works 30[ 3l| where 
the lattice spin configuration is frozen while calculating 
the resistivity, we use here several thousands of configura- 
tions in each of which the resistivity is averaged with many 
thousands of passages. In our previous works, though the 
overall number of MC steps per spin was as high as in 
this work, the fact that we have used only a dozen lattice 
configurations for resistivity calculation has shown strong 
fluctuations in the result. In this work, we have made a 
new device in two steps: 

i) For each lattice configuration all itinerant spins move 
through the system during typically one thousand MC 
steps. Then we thermalize again the lattice for several 
thousands of MC steps before continuing the averaging of 
the resistivity for another thousand MC steps per spin. We 
repeat this cycle for 200 times. In doing so, each itinerant 
spin was averaged over 10^ MC step using " uncorrelated" 
200 lattice configurations in all. 

ii) The resistivity R is defined as = ^ where n is 
the number of itinerant spins crossing a unit area perpen- 
dicular to the X direction per unit of MC time. To know 
this number, we count them at three "detector" surfaces 
perpendicular to the x direction: the first at N^/^^ the 
second at Nx/2 and the third at 3A^a;/4. Averaging the re- 
sistivity at these three system positions helps to improve 
further the results (in our previous works [31| we counted 
them only at the end of the sample). 

As will be shown below, these extensive configuration 
and space averages give much better results with respect to 
those in previous works. The dynamics of itinerant spins 
is governed by the following interactions: 

i) an electric field E which is applied in the x direction. 
Its energy is given by 



-eE- 



(6) 



where is the distance traveled by the itinerant spin 
in a MC step, e its charge. The orientation of is taken 
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at random, its magnitude is taken from a uniform distri- 
bution between and tq where tq is the nearest-neighbor 
distance; 

ii) a chemical potential term (" concentration gradient" 
effect) given by 

nc = Dn{r), (7) 

where n(r) is the concentration of itinerant spins in a 
sphere of radius D2 centered at r. D is a constant taken 
equal to Kq for simplicity; 

iii) interactions between a given itinerant spin and lat- 
tice spins inside a sphere of radius Di (Eq. [3]); 

iv) interactions between a given itinerant spin and other 
itinerant spins inside a sphere of radius D2 (Eq. [2j). 

The simulation is carried out as follows: at a given T 
we calculate the energy of an itinerant spin by taking 
into account all the interactions described above. Then 
we tentatively move the spin under consideration to a new 
position with a step in an arbitrary direction. Note that 
this move is immediately rejected if the new position is in- 
side a sphere of radius Aq centered at a lattice spin or an 
itinerant spin. This excluded space represents the Pauli 
exclusion principle in the one hand, and the interaction 
with lattice phonons on the other hand. If the new posi- 
tion does not lie in a forbidden region of space, then the 
move is accepted with a probability given by the standard 
Metropolis algorithm[32i|. 




Figure 1: Lattice magnetization versus temperature T ior Nz = 8, 
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Figure 2: Resistivity R in arbitrary unit versus temperature T for 
several values of /q: 2 (black circles), 1 (void circles), 0.5 (black 
triangles). Other parameters: Nx = Ny = 20, Nz=8, E = 1, Kq = 
0.5, D = 0.5, Di = 1. 



3. Results on Ferromagnetic Thin Films 

We let A^o itinerant spins travel through the system sev- 
eral thousands of passages until a steady state is reached 
before averaging the spin resistivity. 

The parameters we use in our calculations for ferro- 
magnets are s = S = 1 and = Ny = 20 and = 8. 
Other parameters are Di = D2 = 1 (all distances are in 
unit of the FCC ceh length), Kq = Iq = 0.5 and D = 0.5 
unless otherwise stated, Nq = S x 20^ (namely one itiner- 
ant spin per FCC unit cell), Aq = 0.05 and ro = a/2/2, the 
FCC nearest-neighbor distance. At each T the equilibra- 
tion time for the lattice spins lies around 10^ MC steps per 
spin and we compute statistical averages for the resistiv- 
ity over 10^ MC steps per spin. Note that, as described in 
subsection 2.2, in this work the averaging length 10^ MC 
steps per spin has been divided into 200 segments of 5000 
MC steps; between two consecutive segments we thermal- 
ize again our lattice over several thousands of MC steps 
par spin to explore a maximum of lattice configurations 
encountered by itinerant spins. In doing so we reduce sta- 
tistical fiuctuations observed in our old works |30l l3lj. 

We show in Fig. [1] the lattice magnetization versus T 
for = S, N:, = Ny = 20, J = 1. We find - 9.58 for 
the critical temperature of the lattice spins. 

We show in Fig. [2] the resistivity versus T without 
magnetic field for several values of Iq , interaction between 
itinerant spins and lattice spins. Several remarks are in 
order: 



i) As seen here, no significant fiuctuations of the data 
are observed in the whole temperature range, thanks to 
our new averaging device (see subsection 2.2). 

ii) At Tc, R exhibits a peak at the transition temper- 
ature. The height of the peak decreases with decreasing 
Iq. We see thus that the peak is a consequence of the in- 
teraction between itinerant spins and lattice spins. The 
resistivity shows almost no peak for /q = 0.5. This case 
corresponds to a metal where the interaction between itin- 
erant and lattice spins is very weak. 

iii) We can explain the existence of the peak by the fol- 
lowing argument: the peak is due to the coupling through 
Iq of itinerant spins to the fiuctuations of the lattice spins 
in the critical region around Tc. In our recent work, we 
found from our MC simulation jsij that the resistivity's 
peak is due to the scattering by antiparallel-spin clusters 
which exist when one enters the critical region. Below the 
transition temperature, there exists a single large cluster 
of lattice spins with some isolated "defects" (i. e. clusters 
of antiparallel spins), so that the resistance decreases with 
decreasing T just after Tc. 

iv) However, at very low T, the resistivity increases 
with decreasing T. The origin of this behavior comes from 
the freezing of the itinerant spins due to their interaction 
with lattice spins and with themselves. This is very similar 
to the crystallization of interacting particles at low T. We 
have tested this interpretation by reducing the strength of 
the interactions Kq and Iq. As a matter of fact, R increases 
more slowly with decreasing T. This is seen in Fig. [2]at low 
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Figure 3: Travel path of an itinerant spin at T = 5 (upper) and at 
T = 9.79 (lower) for Iq = 2. Other parameters are the same as in 
Fig. [2] See text for comments. 



T where R is smaller for smaller /q . Note that the increase 
of R at very low T was observed in many experiments on 
various materials not limited to ferromagnets[10-12, 14]. 

v) In the paramagnetic phase, as T increases, small 
clusters will be broken more and more into single disor- 
dered spins, so that there is no more energy barrier be- 
tween successive positions of itinerant spins on their tra- 
jectory. The resistance, though high, is thus decreasing 
with increasing T and saturated as T ^ 00. 

During the simulation, we have followed one itinerant 
spin among 3200 and recorded its successive positions. We 
show in Fig. [3] its travel path at T = 5 where R is very 
low and at T = 9.79 where R is highest. As seen, even 
at T = 5, the spin spends a lot of time to overcome the 
scattering by lattice spins. Almost four thousands of trial 
moves are needed to get through the system. Meanwhile, 
at T = 9.79, the spin under consideration spends all this 
time in a small area. Note that at this peak's temperature, 
very few spins can go until the end unless we increase the 
electric field. 

3.1. Effect of magnetic field. 

Kataoka[22] has shown by a Boltzmann's equation for- 
malism that the magnetic field reduces the peak's height. 
This is what we observed in simulations. We show re- 
sults of R for several fields in Fig. HJ The peak reduction 
is stronger for stronger fields. This is easily understood: 
when a magnetic field is applied on a ferromagnet, the 
phase transition is suppressed because the magnetization 
is not zero at any T. The field reduces critical fiuct nations, 
and hence the number of clusters of antiparallel spins. The 
peak of the resistivity is therefore reduced and disappears 
at high fields. 

To close this section, let us emphasize that in ferromag- 
net s, we have improved the previous results [si^ by using a 
new averaging procedure. We found that the height of the 
peak is intimately related to the strength of the interaction 
between itinerant spins and lattice spins (Fig. [2]). 
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Figure 4: Resistivity R in arbitrary unit versus temperature T, for 
different magnetic fields B: (black circles), 0.25 (void circles), 0.5 
(black triangles), 0.75 (void triangles). Iq = 2 and other parameters 
taken the same as in Fig. 2. 



4. Results on Antiferromagnetic Thin Films 

In the case of antiferromagnets, we study first a film 
with simple cubic (SC) lattice structure. This is because 
the FCC lattice used in the ferromagnetic case shown above 
becomes fully frustrated if we use an antiferromagnetic in- 
teraction. The frustrated case is very particular [33], it 
cannot be treated on the same footing as the non frus- 
trated case. 

Before showing the results on a SC antiferromagnet, 
let us emphasize the following point. The picture of de- 
fect clusters of down spins embedded in a up-spin sea that 
we used above to explain the behavior of the resistivity 
in ferromagnets should be modified in the case of antifer- 
romagnets: in antiferromagnets defects are domain walls, 
clusters on the two sides of a wall both have antiferro- 
magnetic ordering with opposite parity. An itinerant spin 
crossing a domain wall does not have the same scattering 
as in a ferromagnet. Its scattering depends on the num- 
bers of up spins and down spins in the sphere of radius Di . 
In other words, the scattering depends on the energy land- 
scape in the crystal: the itinerant spin will stay a longer 
time where its energy is low, and a shorter time where its 
energy is high. 

The resistivity versus T in zero magnetic field is shown 
in Fig. [5] with Di = D2 = 1. Several remarks are in 
order: i) one observes the absence of a peak of R; ii) the 
variation of R with T has the same shape as the internal 
energy versus T shown in Fig. [6l therefore dR/dT shows 
a peak similar to the specific heat. The peak of dR/dT 
has been experimentally observed in many materials, in 
particular in MnSi [3|, |^ among others [2] . 

The absence of a peak at Tn observed here certainly 
comes from the fact that the motion of an itinerant elec- 
tron is not sharply slowed down at Tn by numerous clus- 
ters of opposite spins. Let us say it again in another man- 
ner: the absence of a peak at the transition is due to the 
fact that the motion of an itinerant spin depends on its 
immediate environment: in ferromagnets, the variation of 
its energy going from a "parallel" cluster to a nearby 
"defect" (or antiparallel) cluster is much larger than the 
energy variation going from a cluster of antiferromagnetic 
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Figure 5: SC AF case. Resistivity R in arbitrary unit versus 
temperature T, in zero magnetic field, with electric field E = 1, 
Io = Ko = 0.5. 
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Figure 7: SC AF case. Resistivity R in arbitrary unit versus Di at 
T = 1, in zero magnetic field, with electric field E = 1, Iq = Kq = 
0.5. 




Figure 6: SC AF case. Internal energy S in unit |J| = 1 versus 
temperature T, in zero magnetic field, with electric field E = 1, 
Io = Ko = 0.5. 
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Figure 8: SC AF case. Resistivity R versus T for several values of 
Di: 1 (black circles), 1.2 (void triangles), 1.4 (black triangles), 1.6 
(void circles), E = 1, B = 0, Iq = Kq = 0.5, D = 0.35. 



ordering to a cluster which is a defect but a defect with an 
antiferromagnetic structure in the SC antiferromagnetic 
case. The smaher AE gives rise to a larger spin mobility 
i. e. a smaller R. Note that experimental data show just 
a shoulder in antiferromagnetic LaFeAsO[16]. 

We show now in Fig. [Tithe effect of Di on the resistiv- 
ity at a given temperature. We observe here an oscillatory 
behavior of R. By analyzing the ratio of numbers of up 
spins and down spins in the sphere of radius Di , we found 
that this ratio oscillates with varying Di: the maxima 
(minima) of R correspond to the largest (smallest) num- 
bers of parallel (antiparallel) spins in the sphere. This 
finding is consistent with what we said before, namely R 
is large when the energy of itinerant spin is low (i. e. large 
number of parallel spins). 

We show in Fig. [8] the resistivity versus T for several 
Di. As seen here, the change of at a given T > Tn with 
varying Di is much smaller than that for T < T^. It is 
interesting to note that R does not depend on Di at T^. 
Further analysis should be carried out to understand this 
behavior at the transition point. 

Finally, to compare with the SC AF case, let us show 
the results of the BCC and diamond-lattice antiferromag- 
nets in Figs. [9]and[T0l As seen, these cases shows no peak 
for any value of Di . Note that R in both cases increases as 
T decreases to zero, unlike the SC AF case. The origin of 
this increase lies in the freezing of itinerant spins at low T. 
The degree of freezing depends on the lattice structure and 
on the strength of the interactions of the itinerant spins 



with their environment as discussed in the ferromagnetic 
section. 

4.1. Discussion 

In the case of ferromagnets, the coupling of the motion 
of itinerant spins to the correlation of the lattice spins gives 
rise to the peak of the resistivity in the transition region. 
Depending on the strength of this coupling, the peak can 
be very sharp or rounded at Tc. The picture of scattering 
by clusters suggested by us is naturally consistent with the 
correlation interpretation. 

In the case of antiferromagnets, the polarized itinerant 
spins are coupled to both parallel and antiparallel lattice 
spins. Due to the opposite correlation signs, their respec- 
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Figure 9: BCC AF case. Resistivity R versus T for several values of 
Di: 1 (black circles), 1.2 (void triangles), 1.4 (black triangles), 1.6 
(void circles), E = 1, B = 0, Iq = Kq = 1, D = 0.5. 
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Figure 10: Diamond-lattice AF case. Resistivity R versus T for Di: 
0.5 (black triangles), 0.75 (black circles), 1 (void circles). E = 1, 
B = 0, Io = Ko = 0.5, D = 0.35. 



tive effects are partially canceled out giving rise to an ef- 
fective coupling weaker than that in ferromagnets. It is 
therefore not surprising that the peak is absent unlike in 
the ferromagnetic case. 

5. Conclusion 

In this work, we have improved our MC simulations 
by averaging the resistivity over a large number of lattice 
spin configurations. The results shown above are much 
better than those in our previous works [30|, HH . Though 
the physics in the ferromagnetic case are not qualitatively 
altered but the precision on the peak position of R is ex- 
cellent and the statistical fiuctuations in the paramagnetic 
phase are reduced. The spin resistivity is strongly depen- 
dent on the temperature. In ferromagnets, at very low T 
the itinerant spins are somewhat frozen. As T increases, 
their motion is thermally activated making a decrease of 
R. However, as T increases further, the system enters the 
transition region, R increases and undergoes a huge peak 
at the ferromagnetic transition temperature. At higher 
temperatures, the lattice spins are disordered, the resis- 
tivity is still large but it decreases with increasing T. The 
existence of the peak in ferromagnets is in a gree ment with 
theories, in particular those of Zarand et al [23| and Haas 
[iT] where interaction between itinerant spins and lattice 
spins is dominant. 

We have also shown here results of some antiferromag- 
nets. The absence of a peak of R in the SC, BCC and 
diamond lattices confirms the prediction of Haas . Let 
us emphasize that our results on frustrated FCC AF [33] 
and on Heisenberg BCC AF [34] show that the shape of 
the resistivity in antiferromagnets strongly depends on the 
spin model and the nature, i. e. first or second order, of 
the lattice phase transition. The extension of the Boltz- 
mann's theory to the case of antiferromagnets is under 
way. 
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course of the present work. 
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